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Abstract. We show that for large instances the order of magnitude of 
the expected Frobenius number is (up to a constant depending only on 
the dimension) given by its lower bound. 



1. Introduction 

Let o be a positive integral re-dimensional primitive vector, i.e., a = 
(ai,...,o n ) T £ Z™ with gcd(a) := gcd(oi, . . . , a n ) = 1. The Frobenius 
number of a, denoted by F (a), is the largest number which cannot be rep- 
resented as a non- negative integral combination of the Oj's, i.e., 

F (a) = max{6 eZ:b^(a,z) for all z G Z^ }, 

where (-, ■) denotes the standard inner product on M. n . In other words, 
F (a) is the maximal right hand side b, such that the well-known knapsack 
polytope P(a, b) = {x £ M> : (a, x) = b} does not contain an integral point. 
From that point of view it is also apparent that the Frobenius number plays 
an important role in the analysis of integer programming algorithms (see, 
e.g., [H El EE! EH EI] ) and, vice versa, integer programming algorithms are 
known to be an effective tool for computing the Frobenius number (see, e.g., 
[HI [131 120] ). There is a rich literature on Frobenius numbers, and for an 
impressive survey on the history and the different aspects of the problem we 
refer to the book [2]. 

Here we just want to mention that only for n = 2 an explicit formula is 
known, which was proven by Curran Sharp in 1884 answering a question by 
Sylvester: 

F (a) = ai a 2 - (ai + a 2 ). 

There is a huge variety of upper bounds on F (a). They all share the prop- 
erty that in the worst case they are of quadratic order with respect to the 
maximum norm of a, say, which will be denoted by |a|oo- For instance, 
assuming a\ < a 2 < . . . < a n , a classical upper bound due to Erdos and 
Graham [14] says 



F(a) < 2a n 
i 
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and, in a recent paper, Fukshansky&Robins [151 (29)] gave an upper bound 
which is also symmetric in the a^s 



(1.1) F(o)< 



(n-l)7r(f + 1) 



W 2 

i=l 



where | • I2 denotes the Euclidean norm. The worst case in all the known 
upper bounds is achieved when the a^s are approximately of the same size, 
and it is also known that in these cases the quadratic order of an upper 
bound cannot be lowered (see HH ;22j). 

On the other hand, Aliev&Gruber [3] recently found an optimal lower 
bound for the Frobenius number which implies that 

(1.2) F (a) > {n - l)!==i (01 a 2 ■ . . . • a n )^ - (ai H h a n ) . 

Hence, if all the aj's are of the same size then the lower bound is only of order 
|a|oo~ ^ ■ In fact, taking the quotient of the (symmetric) upper bound 
(jl.ip with (jl.2p . we see that there is always a gap of order |a|L ■ 

Thus the next natural and important question is to get information on 
the Frobenius number of a "typical" vector o. This problem appears to be 
hard, and to the best of our knowledge it has firstly been systematically 
investigated by V. I. Arnold, see, e.g., [SJ [Tl E] • In particular, he conjectured 
that F (a) grows like j ,1 + 1 /( n - 1 ) for a "typical" vector a with 1-norm |a|i = T 
[5], and in [6j 2003-5] he conjectures that the "average behavior" is 



(1.3) F (a) ~ (n - 1)!^ (ai a 2 • ■ ■ ■ ■ a„)« 



i.e., it is essentially the lower bound. A similar conjecture for the 3-dimen- 
sional case was proposed by Davison [12] and recently proved by Shur, Sinai 
and Ustinov [23]. Extensive computations support conjecture (jl.3p (see [9]). 

In [TU], Bourgain and Sinai proved a statement in the spirit of these 
conjectures, which says, roughly speaking, that 



Prob (F(a)/T w ^ n -^ > D^j < e(D), 



where Prob(-) is meant with respect to the uniform distribution among all 
points in the set 

G{T) = {a e : gcd(a) = 1, |a|oo < T}. 

Here the number e(D) does not depend on T and tends to zero as D ap- 
proaches infinity. The paper [4] gives more precise information about the 
order of decay of the function e(D). Their main result [3J Theorem 1.1] 
implies that 



(1.4) Prob (F (a)/^ 1 ^- 1 ' > D) « n D~ 2 , 

where <C n denotes the Vinogradov symbol with the constant depending on 
n only. In particular, from that result the authors get a statement about 
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the average Frobenius number, namely [U Corollary 1.2] 

T *W) <<nL 

So the order of the average Frobenius number is (essentially) not bigger than 
|a|oc^ 1/ ^ n , which is close to the sharp lower bound (II. 2|) . but there is still 



a gap. 

The main purpose of this note is to fill that gap. We will show that 
Theorem 1. Let n > 3. Then 

Prob fF(o)/(oio 2 • ••• ■ a n )^ > d) <C„ D~ 2 ^ . 



From this result we will derive the desired statement 
Corollary 1. Let n > 3. Then 



EaeG(T) F («)/(«! a 2 ■ 



a 



n) 



1 

I n-1 



T " # G(T) — «»" 1 

These statements fit also perfectly to recent results on the limit distribu- 
tion of Frobenius numbers due to Shur, Sinai and Ustinov [23] and Marklof 
|19j . For instance, in our special setting, |19^ Theorem 1] says that 



lim Prob F (a)/ {a x a 2 ■ . . . • a n )^ > D) = *(£>), 

where : M>o — > M>o is a non- increasing function with ^(0) = 1. 

The proof of Theorem[T]is based on a discrete inverse arithmetic-geometric 
mean inequality which might be of some interests in its own. It will be stated 
and proved in Section 2. Finally, Section 3 contains the proofs of Theorem 
[T]and Corollary [TJ 

2. Reverse discrete AGM-inequality 

For x S R™q the Arithmetic-Geometric-Mean (AGM) inequality states 

that (13^=1 x i) l ^ n — h Y17=i x i- -^ * s known that the "reverse" AGM in- 
equality holds with high probability. More precisely, Gluskin&Milman |16j 
have shown that 



n Ya=1 x 



p ro b 5 „-i v n 1 = -r >a \ < c n a- n ' 2 , 

V(nr=ii^i) 1/n vMnr=iM) 1/n )~ 

where c is an absolute constant and Probgn-iO is meant with respect 
to standard rotation invariant measure on the re-dimensional unit sphere 
S n ~ l = {x S M n : Y17=i x \ = H ere we show an analogous statement 
with respect to the primitive lattice points in the set G(T). 
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With respect to the uniform distribution on G(T), let Lt ■ G(T) — > M>o 
be the random variable defined by 



(n?=i a i 



1/n " 



Theorem 2. Let a > 1 and let k £ {1, ••• ,n — 1}. Then there exists a 
constant c(k, n) depending only on k, n, such that 

Prob (L T >a)< c(k, n) a~ k . 



By Markov's inequality the theorem is an immediate consequence of the 
next statement about the expectation E(L^) of higher moments of Lt- 

Lemma 1. Let k £ {l,...,n— 1}. Then there exists a constant c(k,n) 
depending on k, n, such that 

E(L T ) < c(k,n). 



Proof. First we note that there is an absolute constant c such that 
(2.1) #G(T) > cT n . 

This follows easily from well-known relations between integration and count- 
ing primitive lattice points (see e.g. [HI p. 183, (1)]), but in order to keep 
the paper self-contained we give a short argument here: Let n > 2 and let 

G 2 {T) = {{a, b) T e Z|j : 1 < a, b < T, gcd(a, b) = 1}. 

Since G*2(T) x {1, ... ,T} n ~ 2 C G(T) it suffices to prove the statement for 
n = 2, i.e., G*2(T). There are at most (T/m) 2 pairs (a, b) T G {1, . . . , T} with 
gcd(a, b) = m. Thus 

#G 2 (T) >T 2 (l - fl m j = (2-^ 2 /6)T 2 , 

V m=2 J 

which gives (|2.1D . 
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Now we have 
(2.2) 

1 i / l V^™ 

«0*> = =5r? E = E 



^ 1 ' a€G(T) W K ' aeG(T) V 11 li=l 



1^1 ^ ( k \a\ 1 a l 2 2 



#G{T) ^ n fc 2 \ /,. /'..••• . /„/ {Tr » „ 

^ V aGG(T) j 1 +i 2 +...+i n =fc v z ' ' n/ llli=l a i 

- 1 e if,,,*. ,. ) n e <-''- l/ " 



#G(T) f-^ n k \i 1 ,i 2 ,--- ,i n J H 3 

H+«2H h«n = fc J = l a J =1 

Since for k < n the sum E^™ = i m~ k / n is bounded from above by c(k, n) T 1 ~ k / n , 
where c(/c, n) is a constant depending only on k and n, we find 

T 

(2.3) Yl a ?~ fc/n ^ c ( fc > ^)^ 1+ij ~ fc/n - 

aj=l 

Thus, for «i + «2 H r- % = k we obtain 

n^ay- fc/n <c(fc,n) ft T» 

j = l CLj=l 

Hence we can continue (|2.2p by 

E (4)<^ E >( jl ,J..,J*" , " T " 

*1+*2H Mra = fc 

Tin 

Finally, with (|2.ip we get the assertion. □ 

Remark 1. With a little more work one can prove Lemma[J\for any positive 
number k < n, not only for integers. 

Next we want to point out that the arguments in the proof of Lemma Q] 
lead to the following lower bound on the random variable X? '■ G(T) — > R>o 
given by 

F(a) 



X T (a) 



i_ 



Proposition 1. Let n > 2. T/ien i/iere exists a constant depending only on 
n such that 

E(X T ) > (n - 1)!^ (l - c(n) 
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Proof. On account of the lower bound (II. 2p on F (a) it remains to show that 



i £ F =1 \_ < c (n) r-^i. 



Following the argumentation in (|2.2p we find 

E j^r ^ E E II 



aeG(T) (nj=l a «)" _1 U+12H Nn=l ai,a 2 ,...,a n =l j=l 

T 

n 

ai,a2,...,a n =l 



ai n " 1 a 2 n - 1 •...•a ri "- 1 



r ^ , r ^ n-1 

n > a, V a. ' 1 1 



ai=l \a2=l / 

Thus, analogously to (|2.3j) . and on account of (12. ip we obtain 

# a <r> J? m <xiu^ # G < T » 1 

1 

< c(n)T n ~ 1 . 



□ 



3. Proof of Theorem [T] and Corollary Q] 

We keep the notation of the previous section. First we note that (jl.4p is 
certainly also true for any other norm in the denominator, in particular for 
the 1-norm | • |i- Thus 

(3.1) Prob [JLg- >£>) <£„ / 



lab - 1 



Secondly, we observe that 
(3.2) 



Prob | w"- 1, , > J = Prob ( aw')-, ji) h t 



(ai • . . . • a™)™- 1 / • • • • • (»ti)"- 



■n 



Pro], | L^" 1 > ( - " 



1 „_A (n-1) 2 

Prob [ Lt > —7 n I <c n 7 
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where in the last inequality we have applied Theorem [2] with k = n — 1. 
Together with (|3.1j) we obtain 



Prob (X T (a) > (3) = Prob 



i 

F(a) -m ' 




for any t G (0, 1). With t = (n — l)/(n+ 1) we finally get 

o n — 1 

Prob (X T (a) >/5) « n p~ 2 ~ , 

which shows Theorem [TJ 

For the proof of Corollary [I] we note that 



E(X 7 



roo roc 

J Prob (X T > x) dx < 1 + y Prob (X T > x) dx 



poo , 

o n — 1 

<C n 1 + / x "+i dx. 



For n > 4, the last integral is finite and so we have 
(3.3) E(X T ) « n 1. 

For the case n = 3 we just note that on account of Remark [Done can also 
bound Prob (Xr(a) > (3) by a function like (3~ <yl+t ' > and so we also get (|3.3p 
in this case. Together with Proposition [Q Corollary [1] is proven. 
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